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Abstract In various astrophysical contexts, we analyze 
self-similar behaviours of magnetohydrodynamic (MHD) 
evolution of a quasi-spherical polytropic magnetized gas 
under self-gravity with the specific entropy conserved 
along streamlines. In particular, this MHD model analy- 
sis frees the scaling parameter n in the conventional poly- 
tropic self-similar transformation from the constraint of 
n + 7 = 2 with 7 being the polytropic index and there- 
fore substantially generalizes earlier analysis results on 
polytropic gas dynamics that has a constant specific en- 
tropy everywhere in space at all time. On the basis of the 
self-similar nonlinear MHD ordinary differential equa- 
tions, we examine behaviours of the magnetosonic criti- 
cal curves, the MHD shock conditions, and various asymp- 
totic solutions. We then construct global semi-complete 
self-similar MHD solutions using a combination of ana- 
lytical and numerical means and indicate plausible astro- 
physical applications of these magnetized flow solutions 
with or without MHD shocks. 
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1 Introduction 

Dynamic evolutions of various astrophysical shock flows 
under gravity on completely different spatial and tem- 
poral scales have been attracting astrophysicists, mainly 
because of a rich variety of observational phenomenol- 
ogy in different astrophysical systems, such as stellar 
collapses, supernova explosions and supernova remnants 
(e.g., Bethe et al. 1979; Goldreich & Weber 1980; Cheva- 
lier 1982; Yahil 1983; Lattimer & Prakash 2004; Lou & 
Wang 2006, 2007; Lou & Cao 2008; Lou & Hu 2008 in 
preparation), star formation, proto-stellar core forma- 
tion as well as 'champagne flows' in star-forming molec- 
ular clouds (e.g., Shu 1977; Shu et al. 1987; Tsai & Shu 
1995; Shu et al. 2002; Shen & Lou 2004; Bian & Lou 
2005; Lou & Gao 2006), formation and spherical accre- 
tions of black holes (e.g., Bondi 1952; Bahcall & Ostriker 
1975; Hennawi & Ostriker 2002; Cai & Shu 2005; Hu et 
al. 2006), active galactic nuclei (AGNs) at the galactic 
level (e.g., Small & Blandford 1992) and the formation of 
galaxy clusters and voids in the Universe (e.g., Gunn & 
Gott 1972; Fillmore & Goldreich 1984a, b; Bertschinger 
1985) as well as dynamic galaxy cluster winds (Lou, 
Jiang & Jin 2008; Jiang & Lou 2008 in preparation). 
This broad class of theoretical models poses challenges 
for a deeper and simpler theoretical understanding of 
relevant astrophysical shock flows under consideration. 
Meanwhile, these problems also involve a diverse range of 
physical processes and the relevant model investigations 
often involve numerical hydrodynamic simulations with 
data input from high-energy physics experiments and/or 
theories (e.g., Janka & Miiller 1996 for supernovae). 

In this paper, we mainly focus on basic hydrodynamic 
and magnetohydrodynamic (MHD) aspects of this class 
of problems because of the relative simplicity and gener- 
ality, allowing for various novel and interesting solution 
features to manifest. We introduce the general polytropic 
description to subsume several unspecified energetic pro- 
cesses into a few index parameters. There might be cer- 
tain evolution phases where a simple and direct anal- 



ysis works almost as well as numerical simulations do 
(e.g., Janka & Miiller 1996), and in this situation such 
an analysis would be extremely valuable for physical in- 
sight (e.g., Yahil 1983) and for testing numerical codes 
in construction. In this spirit, we advance a theoretical 
MHD flow model with this kind of simplification (e.g., 
self-similarity) and idealization (e.g., a completely ran- 
dom transverse magnetic field and specific entropy con- 
servation along streamlines) to reveal plausible MHD be- 
haviours and shocks of magnetized gas flows of astro- 
physical interests. 

While self-similar evolutions of gas dynamics in spher- 
ical geometry have been actively pursued early on (e.g., 
Sedov 1959; Bodcnhcimcr & Swcigart 1968; Larson 1969a, 
b; Penston 1969a, b), extensively analyzed by many in 
diverse astrophysical contexts (e.g., Shu 1977; Hunter 
1977; Cheng 1978; Goldreich & Weber 1980; Chevalier 
1982; Yahil 1983; Whitworth & Summers 1985; Tsai & 
Hsu 1995; Shu et al. 2002; Shcn & Lou 2004; Fatuzzo ct 
al. 2004; Lou & Gao 2006) and significantly reformulated 
in various physical aspects (e.g., Terebey, Shu & Cassen 
1984; Suto & Silk 1988; Chiueh & Chou 1994; McLaugh- 
lin & Putritz 1997; Boily & Lynden-Bell 1995; Cai & Shu 
2005), we have found novel features and generalized sim- 
ilarity solutions for this type of classical nonlinear prob- 
lems in recent years. As examples, our new similarity 
solution features include envelope expansion with core 
collapse (EECC) solutions which smoothly pass through 
the sonic critical line twice (Lou & Shen 2004; Lou & Gao 
2006), quasi-static asymptotic solutions (Lou & Wang 
2006), quasi-static asymptotic MHD solutions (Lou & 
Wang 2007; Wang & Lou 2007), quasi-static asymptotic 
solutions in two gravity-coupled fluids (Lou, Jiang & Jin 
2008; Jiang & Lou 2008 in preparation), strong magnetic 
field solutions (Yu & Lou 2006; Wang & Lou 2007), and 
various new (MHD) shock solutions (Shcn & Lou 2004; 
Bian & Lou 2005; Yu ct al. 2007). 

In addition to the known Larson-Penston (LP) type 
solutions and central free-fall solution (Shu 1977), the ex- 
istence of several new nonlinear self-similar MHD poly- 
tropic flow solutions reveals multiple clues and leads to 
new concepts. 

First of all, knowing these self-similar nonlinear solu- 
tions, we realize that a number of self-similar behaviours 
may possibly occur in seemingly similar flow systems, 
and the evolution problem of 'which goes where' remains 
an important open question. Clearly, this cannot be im- 
mediately answered within the self-similar framework. 
While it is generally thought that self-similar behaviour 
gradually emerges when a dynamical system evolves suf- 
ficiently far away from its initial and boundary condi- 
tions, we would emphasize that initial and boundary con- 
ditions do affect the nonlinear self-similar behaviour in 
the sense of which similarity behaviour emerges even- 
tually. In other words, the final self-similar phase still 
retains memory of a certain class of initial and bound- 
ary conditions. As an example, Tsai & Hsu (1995) intro- 



duced a self-similar behaviour that has the outer (initial) 
configuration the same as that of Shu (1977), but their 
inner (final) behaviour is qualitatively different from that 
of Shu (1977): the model of Tsai & Hsu (1995) includes 
a self-similar isothermal shock and also the central mass 
accretion rate differs from that of Shu (1977). As sug- 
gested by simulation results of Tsai & Hsu (1995), a 
stronger push at the centre of a static isothermal sphere 
may result in a shock, while a weaker push will lead to 
an isothermal expansion-wave collapse solution (EWCS; 
see their Figs. 1 — 3). 

Secondly, although former nonlinear self-similar so- 
lutions were often obtained from numerical simulations 
first (e.g., Larson 1969a, b; Penston 1969a, b), the very 
existence of different self-similar behaviours in the same 
physical setting would imply that so far numerical hy- 
drodynamic simulations might have not been sufficiently 
thorough in exploring the overall dynamic evolution fea- 
tures of nonlinear flows with or without shocks. These 
semi-analytical solutions are important clues and bench- 
marks for numerical code development. More importantly, 
numerical simulations can further tell whether new self- 
similarity solutions obtained in our analysis can be actu- 
ally realized under realistic flow situations, by introduc- 
ing different initial and boundary conditions for a similar 
flow system and by observing which similarity behaviour 
does it manifest eventually. 

Finally, it is also possible that one flow system may 
not behave self-similarly in the strict sense, nonetheless 
its evolution can be qualitatively explained by intuitions 
gained by examining various self-similar solutions: for ex- 
ample, the central free-fall phase of Shu (1977) is com- 
mon for a wide range of polytropic index 7; the existence 
of other new solutions requires specialized conditions. 
Physically, our new asymptotic solutions have different 
features in terms of the balance and competition of the 
involved forces, which tells us about the possible asymp- 
totic flow behaviours. 

As an essential ingredient in our theoretical model 
development, magnetic fields are important in various 
astrophysical settings; their presence in ionized plasma 
systems is ubiquitous and they play dynamical as well as 
diagnostic roles in many ways; they are responsible for 
sporadic violent activities and for producing relativistic 
particles such as cosmic rays through shocks and recon- 
nection. The generation of magnetic fields is related to 
convective motions, sustained turbulence and differential 
rotations, via dynamo effects or magnetorotational insta- 
bilities (MRI) (e.g., Parker 1979; Thompson & Duncan 
1993 for dynamo effects of magnetized pulsars; Chan- 
drasekhar 1961; Balbus & Hawley 1998; Balbus 2003 for 
disc magnetorotational instabilities). For the purpose of 
our formulation here, we take the presence of magnetic 
field for granted and simply presume that a completely 
random magnetic field permeates in a gas medium with 
self-gravity This gas medium can be stellar interior, hot 



stellar coronae, interstellar medium (ISM), or intraclus- 
ter medium (ICM) and so forth. 

During the collapse phase of a quasi-spherical sys- 
tem under self-gravity, a highly conducting magnetofluid 
undergoes a magnetic field enhancement through the 
magnetic flux conservation, commonly referred to as the 
frozen-in condition for magnetic field. When rotation is 
sufficiently slow in an astrophysical system, the over- 
all geometry may remain quasi-spherical and our model 
analysis with the quasi-spherical random-field approxi- 
mation (e.g., Zel'dovich & Novikov 1971) can be appli- 
cable (Yu & Lou 2005; Yu et al. 2006; Lou & Wang 2007; 
Wang & Lou 2007; Lou & Hu 2008, in preparation). Chi- 
ueh & Chou (1994) first studied a quasi-spherical isother- 
mal MHD problem by including the magnetic pressure 
gradient into the radial momentum equation. However, 
the magnetic tension force does not average to zero (Yu 
& Lou 2005; Wang & Lou 2007) and thus should be 
also included. Although on small scales, magnetic ten- 
sion force tends to drive inhomogeneous gas motions, we 
presume that these non-spherical flows may be neglected 
as compared to the large-scale radial bulk motion of gas. 
The key point is that the (self-) gravity is strong enough 
to hold on the entire gas mass and induce core collapse. 
Hence on large scales in a quasi-spherical geometry, a 
completely random magnetic field contributes to the dy- 
namics in the form of the average magnetic pressure gra- 
dient force and the average magnetic tension force in the 
radial direction. 

In a series of papers during past several years, we 
have systematically investigated isothermal MHD prob- 
lems (Yu & Lou 2005; Yu et al. 2006) and polytropic 
MHD problems in the special case of a constant specific 
entropy in space at all time (i.e., p — np 1 with a constant 
K, a pressure p, a mass density p, and a polytropic in- 
dex 7; Wang & Lou 2007; Lou & Wang 2007; Lou & Hu 
2008, in preparation). Under certain circumstances, the 
isothermality may be a crude yet sensible first approxi- 
mation for astrophysical fluid dynamics with or without 
magnetic field (e.g., for star formation processes within 
a collapsing molecular cloud, Shu 1977; hot bubbles and 
superbubbles in ISM produced by supernovae, Lou & 
Zhai 2008, in preparation; or central core collapse region 
of a globular cluster, Bahcall & Ostriker 1975; Inagaki & 
Lynden-Bell 1983), even for modeling shocks in dynam- 
ical processes (e.g., Courant & Friedrichs 1976; Spitzer 
1978). In other situations, a polytropic fluid with a con- 
stant specific entropy is invoked for theoretical model 
investigation (e.g., Goldreich & Weber 1980 for a stellar 
collapse of a relativistically hot gas priori to the rebound 
process; Yahil 1983 for gravitational core collapses; Suto 
& Silk 1988; Lou & Cao 2008 for rebound processes in 
stellar core collpases). While the polytropic approxima- 
tion with a constant specific entropy is a useful descrip- 
tion for various astrophysical processes, there is however 
no obvious reason why the entropy should necessarily 
remain constant in a polytropic gas. In dynamic pro- 



cesses, a more general situation involves the 'specific en- 
tropy' conservation along streamlines; and this includes 
the constant entropy as a special case. Therefore, self- 
similar polytropic gas evolution under this 'equation of 
state' represents a significant generalization of the 'con- 
stant entropy' problem, and can be more adaptable to 
various dynamic model applications in astrophysics. 

Indeed, several authors have adopted such an 'equa- 
tion of state' and carried out their model analyses to var- 
ious extents (e.g., Cheng 1978; Chevalier 1982; Fatuzzo 
et al. 2004; Wang & Lou 2007; Lou & Cao 2008; Lou 
& Hu 2008, in preparation). In technical terms, adopt- 
ing such an 'equation of state' is equivalent to freeing 
the choice of index parameter n in the usual self-similar 
transformation equation (|7|) from the constraint 71+7 = 2 
with 7 being the polytropic index (Wang & Lou 2007). 
Here, this index parameter n holds the key to the scaling 
of physical variables under consideration, such as mass 
density, random magnetic field strength, radial flow ve- 
locity, pressure, and temperature in the outer portion 
of a quasi-spherical gas. As 7 is a thermodynamic pa- 
rameter, we expect that n and 7 parameters should be 
independent of each other in general. For instance, Lou & 
Wang (2007) have shown that the radial scaling of a com- 
pletely random magnetic field is actually independent of 
the polytropic index 7. The main motivation of this pa- 
per is to revisit the self-similar MHD problem (Wang & 
Lou 2007) but with an equation of state for specific en- 
tropy conservation along streamlines and to generalize 
our new self-similar hydrodynamic and MHD solutions 
into a substantially larger parameter regime. 

This paper with special emphasis on theoretical as- 
pects is structured as follows. Our motivation and back- 
ground information are provided in Section 1 as an in- 
troduction. Section 2 describes the basic formalism on 
nonlinear MHD flows with a completely random trans- 
verse magnetic field, including the self-similar transfor- 
mation and MHD shock conditions. Section 3 presents 
various analytical asymptotic MHD solutions for large 
and small x as well as along the magnetosonic critical 
curve; these asymptotic MHD solutions are valuable for 
understanding relevant physics and can be utilized to 
construct global semi-complete solutions which are valid 
in the range of < x < +00. Section 4 shows examples 
of various possible semi-complete global solutions with 
or without MHD shocks and outlines astrophysical ap- 
plications of these solution types. We discuss our results 
in Section 5. Mathematical details are included in several 
appendices for the convenience of reference. 



2 Formulation of the Model Problem 

2.1 Nonlinear Magnetohydrodynamic Equations 

In spherical polar coordinates (r, 9, 0), the basic non- 
linear equations for a quasi-spherical MHD evolution of 



a general polytropic gas under self-gravity include the 
mass conservation equation 
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where p(r,t) is the gas mass density, u(r,t) is the bulk 
radial flow speed, M(r,t) is the enclosed mass within 
radius r at time i; the radial momentum equation 
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where p is the thermal gas pressure, < B 2 > is the en- 
semble average of the random transverse magnetic field 
B t squared, G — 6.67 x 10" 8 dyne cm 2 g~ 2 is the grav- 
itational constant; the magnetic induction equation 
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and the conservation equation for 'specific entropyLU 
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along streamlines. Note that the Poisson equation relat- 
ing the mass density and the gravitational potential is 
consistently satisfied under the approximation of a quasi- 
spherical symmetry. This set of MHD equations is the 
same as that of Wang & Lou (2007) and Lou & Wang 
(2007), except for equation ©, which is the key differ- 
ence between this analysis and our former analyses. In 
other words, the equation of state p = up 1 with k being 
a constant is only one special case satisfying equation ([6]) 
and is referred to as the 'usual' or 'conventional' poly- 
tropic equation of state. Equation of state ([6|) has been 
adopted in various previous works as noted in Section 1. 
The assumptions on random magnetic field distribution 
(Zel'dovich & Novikov 1971) and magnetic induction are 
the same as those in Yu & Lou (2005), Yu et al. (2005), 
Lou & Wang (2007) and Wang & Lou (2007). The work 
of Yu & Lou (2005) studies the free-fall collapse of an 
isothermal magnetized gas cloud (Shu 1977; Chiueh & 
Chou 1994) and is relevant to star formation and to for- 
mation of hot bubbles as well as superbubbles in ISM. In 

The concept of 'specific entropy' is here extended to sit- 
uations where 7, simply regarded as an index parameter, is 
not necessarily the ratio of actual gas specific heats. 



contexts of shocks and 'champagne flows' (Tsai & Hsu 
1995; Shu et al. 2002; Bian & Lou 2005), the MHD model 
of Yu et al. (2005) extends the scenario and analysis 
to magnetized isothermal clouds. The models of Lou & 
Wang (2006, 2007) are designed for rebound shock pro- 
cess in supernovae without or with a random magnetic 
field and explore the origin of intense magnetic fields 
on compact objects. Given the above approximations, a 
proper combination of equations (TTJ — ^ leads to an 
MHD energy conservation [see equation (7) of Wang & 
Lou (2007) and equation (4.2) of Fan & Lou (1999)]. 

By performing the time reversal transformation t — > 
—t, r — > r, u —* —u p—*p,p—* , p,M^fM and 
< B\ >^>< Bf >, the above MHD equations are invari- 
ant. Using this time reversal invariant property one can 
construct complete global solutions from semi-complete 
global solutions (e.g., Lou & Shen 2004). 



2.2 MHD Self-Similarity Transformation 

In order to transform nonlinear MHD partial differential 
equations (QJ — ([6]) to a set of nonlinear MHD ordinary 
differential equations (ODEs), we introduce the following 
MHD self-similarity transformation, namely 
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Here, a; is a dimensionless independent variable, and v(x), 
a(x), fi{x), m{x) and w(x) are functions of x only; k is a 
dimensional scaling factor to consistently make x, v, a, (3, 
m and w dimensionless; n is an important scaling param- 
eter as noted in Section 1, which essentially controls the 
scaling of dimensional physical quantities with respect 
to time t and radius r. As usual, we refer to v(x), a(x), 
f3(x), m(x) and w(x) as the reduced radial flow speed, 
mass density, thermal gas pressure, enclosed mass, and 
magnetic energy density (associated with the completely 
random transverse magnetic field), respectively. 

Under self-similarity transformation ([7]), equations ([2]) 
and ^ together give the reduced enclosed mass m(x) as 



m = ax 2 (nx — v) , 



(8) 



where nx — v > is required in order to ensure that 
the total enclosed mass remains positive (see Appendix 
[S]for more details). Also under the same transformation 
equation [[?])■ mass conservation equation ([1]) becomes 

(nx — v)a — av — —2(x — v)a/x , (9) 

while magnetic induction equation (O reads 

(nx — v)w' — 2wv = 2[v — (2 — n)x]w/x . (10) 



Equations J9J) and (|10|) together lead to the first integral 
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where the integration constant h is a dimensionless pa- 
rameter representing the strength of random transverse 
magnetic field B t . In essence, relation (fTTj) physically cor- 
responds to the magnetic frozen-in condition. The above 
treatment parallels those of Wang & Lou (2007) and of 
Yu & Lou (2005), and the results are the same. Com- 
bining equations (JSj) and ([9]), the 'specific entropy' con- 
servation equation ([6]) along streamlines is equivalent to 
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where exponent parameter q = 2(n + 7 — 2)/(3n — 2) 
and hence 7 = 2 — n + (3n — 2)q/2. In general, equa- 
tions ([2]) and (HJ) imply that "PP 1 (directly related to 
the 'specific entropy') can be a fairly arbitrary function 
of the enclosed mass M(r, t) along streamlines. It is our 
requirement of self-similarity solutions that leads to the 
specific power-law form of relation (fl2|) . The key point 
of this analysis is that q ^ or n + 7 ^ 2 is allowed 
in general. As noted by Lou & Cao (2008), here we do 
not need a proportional coefficient in equation (|12| for 
7 ^ 4/3. The reason is simply that if we write 



/3 = Ca l m q 
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where C is an arbitrary scaling coefficient, then a change 
of parameter k in self-similarity transformation equation 
(J7J) to C 1/(1 ~ 3?/2) fc would make the coefficient C disap- 
pear. We shall not consider the special case of 7 = 4/3 
or q = 2/3, in which there exists yet another parameter 
as an multiplicative scaling factor between (3 and a y m q 
(see Lou & Cao 2008 for a more detailed analysis on the 
special case of 7 = 4/3 which substantially generalized 
the earlier research analysis of Goldreich & Weber 1980 
and Yahil 1983). Taking equations ©, (QTJ and ([12]) to- 
gether, we obtain the following MHD ODE 

[ a 2-n+3n,/2 I 2, (M _ y yy^ _ ^ _ y y + ^ 2q ,/ 



= -(n - l)v - 2hax - (nx - v)a/(3n - 2) . (14) 

Equations © and (TTJ} together lead to separate equa- 
tions for a' and 1/ respectively in the compact form of 



X(x, a, v)a' = A(x, a, v) , X(x, a, v)v' — V(x, a, v) , 

(15) 

where the three functional coefficients X, A, and V are 
defined explicitly below, namely 
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Equations (JT5]) and (|16p are coupled nonlinear MHD 
ODEs and with specified 'boundary' and 'initial' con- 
ditions (i.e., appropriate asymptotic solutions at large 
and small x) , they can be integrated numerically by the 
standard fourth-order Runge-Kutta method (e.g., Press 
et al. 1986); one needs to pay special attention on the 
magnetosonic singular surface where X(x, a, v) = in 
both ODEs of (fT5|) . The solutions of these two coupled 
nonlinear MHD ODEs cannot smoothly cross the sin- 
gular surface unless they cross it with a MHD shock, 
or they satisfy the critical condition on the so-called 
magnetosonic critical curve where both the denominator 
X(x, a, v) and the numerators A(x,a,v) and V(x, (X, v) 
of these two equations in (|15|) vanish simultaneously. The 
magnetosonic critical curve can be determined numeri- 
cally by the specific procedure described in Appendix 151 
and the magnetosonic critical condition is derived and 
solved numerically in Appendix [C] As necessary checks, 
our analysis here reduces to earlier known results (e.g., 
Shu 1977; Suto & Silk 1988; Lou & Shen 2004; Yu & Lou 
2005; Bian & Lou 2005; Lou & Wang 2006; Wang & Lou 
2007) if we set the relevant parameters to appropriate 
values (e.g., q = 0, 7 = n = 1, and h = 0, and so forth). 



2.3 MHD Shock Conditions in Self-Similarity Form 

An MHD shock is characterized by discontinuities in 
thermal gas pressure, mass density, temperature, tan- 
gential magnetic field and radial flow velocity, caused 
by steepening of nonlinear MHD flow evolution or by 
various explosion and jet processes. In our self-similar 
and quasi-spherical symmetric formulation for a general 
polytropic gas, a self-similar MHD shock can be readily 
constructed mathematically, which represents a possible 
evolution of an MHD shock in such a magnetized flow 
system (see e.g., Chevalier 1982; Tsai &: Hsu 1995; Shu 
et al. 2002; Shen & Lou 2004; Bian & Lou 2004; Yu et al. 
2006; Lou & Wang 2007 for earlier results on self-similar 



shocks either with or without a completely random mag- 
netic field). 

To characterize an MHD shock specifically, we need 
to impose several conservation laws across the two sides 
of a shock front in the comoving framework of reference. 
Denoting u as the flow velocity, u s as the shock velocity, 
p as the gas mass density, p as the thermal gas pressure, 
and < B 2 > as the mean square of the random trans- 
verse magnetic field, respectively, we have in the shock 
framework of reference the mass conservation 
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the MHD energy conservation equation 
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and the magnetic induction equation 
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This set of jump conditions for an MHD shock is the 
same as that adopted in Lou & Wang (2007). We here use 
a pair of square brackets outside each expression enclosed 
to denote the difference between the upstream (marked 
by subscript '1') and downstream (marked by subscript 
'2') quantities, as has been done conventionally for shock 
analyses (Landau & Lifshitz 1959; Zel'dovich & Raizer 
1966, 1967). 

There are two parallel sets of MHD self-similar trans- 
formation in the upstream and downstream domains re- 
spectively, because the parameter k in self-similarity trans- 
formation equation (J7|) is related to the sound speed and 
thus can be different in the two flow regions across a 
shock (see e.g., Lou & Wang 2006 and Lou & Gao 2006). 
We shall set k 2 = \ 2 ki with the scaling ratio A repre- 
senting this difference in the upstream k\ and the down- 
stream k 2 ■ The following relations are necessary for con- 
sistency 
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and we shall simply use h instead of hi and h 2 from 
now on. Using the two relations in (J2TJ) . MHD shock 
jump conditions l|T7|) to (|2"U)) can be readily cast into the 
following self-similarity form 



ai(nxi —Vi) = \a 2 (nx 2 - v 2 ) , 

2-n+3nq/2 In, ^o . /• \2 

a x x x [nx\ — vi) H + a.\{nx\ — v\) 



(22) 



ha 2 x 2 



= y 



2 — n+3nq/2 2o / \o 

i 2 x 2 (nx 2 - V 2 ) H 



+a 2 (nx 2 — v 2 ) 2 + 



«2 X 2 



(23) 
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(24) 



These three self-similar MHD shock equations (|22]l- (|24|) 
can be explicitly solved and relevant details of derivation 
can be found in Appendix [Dl 



3 Asymptotic MHD Solutions 

We connect relevant asymptotic MHD solutions by nu- 
merical integration to construct semi-complete global so- 
lutions for self-similar MHD flows with or without shocks. 
Asymptotic analytical MHD solutions are extremely valu- 
able and illustrate the balance and dominance of various 
forces such as gravity, pressure force and Lorentz force. 
All MHD solutions obtained here are more general with 
q =/= (i.e., n + 7 7^ 2) and are consistent with previ- 
ous results published in the literature (Shu 1977; Suto 
& Silk 1988; Lou & Shen 2004; Yu & Lou 2005; Wang 
& Lou 2007; Lou & Wang 2006, 2007; Lou & Li 2007 in 
preparation; Lou & Cao 2008). 



3.1 MHD Solutions of Pressure Dominance 

When the thermal pressure becomes significant in an 
asymptotic MHD solution, it is generally expected that 
the pressure term will enter the asymptotic ODEs gov- 
erning this solution. Therefore, parameter q will appear 
in the leading order of the asymptotic MHD solution. 

3.1.1 Asymptotic MHD Solutions of Finite Density and 
Velocity at Large x 

For asymptotic MHD solutions at large x, the gravita- 
tional force, magnetic force and thermal pressure gra- 
dient force are all in the same order of magnitude, and 
v(x) and a(x) become finite or approach zero in the limit 
of x — > +00. In this case of large x, we have to leading 
order the following coupled nonlinear MHD ODEs 



a 



2a 

nx 




(n — l)v 

— + 

nx 


1 ( 2h(n- 1)" 
_(3n-2) ' n 2 



a 



-2(2 - n)n"- V-" +3n,z/ V^ 2 



(25) 



and the corresponding asymptotic MHD solution is 



3.1.2 Magnetostatic Solution for a Polytropic Sphere 



Ax- 2 ' n + 



Bx 1 ' 1 ^ 



(26) 



2h(n - 1) 



(3n - 2) 



A 



+2(2-n)n q - 1 A 1 - n+3nq l 2 \x 1 - 2/n + --- , (27) 



where A and B are two constants of integration. The ra- 
dial profile of magnetic field can be readily derived from 
equation (fTTj) . We find that inequality n > 2/3 (i.e., a 
positive enclosed mass M) is sufficient to warrant the va- 
lidity of this asymptotic MHD solution. For 2/3 < n < 1 , 
the two coefficients A and B can be fairly arbitrary. The 
special isothermal case (Shu 1977) corresponds to B = 0, 
q = 0, n = 1, 7 = 1, and h — 0; and the B / 0, (j = 0, 
n=l,7 = l,/i = version was implicit in Hunter 
(1977). The constant coefficient B is a free parameter 
in this solution, because v — Bx 1 ^ 1 / 71 itself would sat- 
isfy du/dt + udu/dr — 0, as if no force is active. For 
n = 1, the flow approaches a constant radial velocity at 
the outer portion (Whitworth & Summers 1985; Lou & 
Shen 2004). This asymptotic solution has been applied 
to various collapse problems (e.g. Shu 1977; Yahil 1983; 
Lou & Shen 2004; Bian & Lou 2005; Lou & Wang 2006; 
Lou & Gao 2006), because of the common feature of 
inflow, outflow, contraction and/or collapse of dynamic 
evolutions. Note that n = 1 case here differs from that 
of Suto & Silk (1988) because we require specific en- 
tropy conservation along streamlines. For 1 < n < 2, we 
should set B = to avoid the flow speed divergence of 
the first term in solution (|2T|) at large x. With n = 2 
and B = 0, we have a constant asymptotic inflow speed 
again by solution (f2"7) . For n > 2, we must require both 
leading order terms in flow speed solution (|27|) to vanish 
to avoid velocity divergence; nevertheless for coefficient 
A > 0, the second term in flow speed solution ([27| does 
not vanish; this second term would vanish for a trivial 
solution of A = 0. Unless other constraints appear, we 
must require 2/3 < n < 2, corresponding to a range of 
mass density scaling p oc r~ 3 to p oc r _1 . For a compar- 
ison, Cheng (1978) indicated the more obvious limit of 
p oc r~ 3 for n — » 2/3, while our analysis here also gives 
another limit of p oc r _1 for n — > 2 due to the second 
term in asymptotic solution (|2"7| for the flow speed. At 
large x, the temperature scales as x 2{jl ~ 1 ^l n . For n < 1, 
the temperature remains finite, while for n > 1, the tem- 
perature diverges as x — > oc and one needs to invoke a 
finite size of an astrophysical system in order to make 
use of this solution. 



Setting v = for all x in equation (|15[) , one immediately 
obtains an exact magnetostatic solutiorQ 



r^2 



2(1 - n)(3n - 2)h 

1 

2(2-n)(3n-2) 



-i -l/(n-3ng/2) 



-2/n 



= A x- 2 / n . (28) 

In the limit of f — > 0, we may still regard equation (jTTj) 
to be valid from the MHD evolutionary perspective; it is 
then straightforward to derive the corresponding radial 
profile of the random magnetic field. This magnetostatic 
solution for a singular polytropic sphere (SPS) is glob- 
ally valid for parameters that physically allow expression 
(f2"5)l except for the central divergence as x — > + . The key 
difference here is that q ^ orn+7 j^ 2 is allowed in gen- 
eral in SPS solution ([2"5]) . The quasi-magnetostatic solu- 
tion described immediately below is constructed based 
on SPS solution (f2"5|) . Of course, all forces present are 
significant in this solution. 

3.1.3 Quasi- Magnetostatic Solutions at Small x 

For quasi-magnetostatic asymptotic solutions (see Lou 
& Wang 2006, 2007 and Wang & Lou 2007 for details of 
deriving the hydrodynamic and MHD counterparts with 
q = 0; see also Lou, Jiang & Jin 2008 and Jiang & Lou 
[2008, in preparation] for this solution in two gravity- 
coupled fluids) at small x values, we first substitute 



A x 



-2/n 



Nx 



K-l-2/n 



Lx 



K 



(29) 



into coupled nonlinear MHD ODEs (TT5")) and (fl7J)) , where 
Aq is explicitly defined by equation (|28p , and then obtain 
two nonlinear algebraic equations for the three relevant 
parameters N , K and L, namely 



n(K - 1)N = (K + 2 - 2/n)A L , 



(30) 
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Here we have introduced a handy parameter 

n 2 (1 — n) 



<l 



2(2-n)(3»-2) (2-n) 



6 )Q \N . (31) 



(32) 



for simplicity and for the convenience of mathematical 
derivation. Equations (|3T)|) and (f3"Tj) together give rise to 

2 In equation (11) of Lou & Wang (2006), there are two 
typos; in the coefficients of both a(x) and m(x) there, the 
exponent should be — 1/n instead of 1/n. 



a quadratic equation for K in terms of three parameters 
n, Q and h, namely 



2(3n-2) 



i ( 3 ™ - 2 ) ^ 
nh + - -Q 



K^^-^K 



+ 



2(2-n)(l 



» 



(3n-2) 2 (l-4/n)Q 
(3n - 2) 



. 
(33) 



Once proper roots of K are known, parameters N and L 
are proportional to each other by equation (|3"0|) and only 
one of them is free to choose. The existence of the exact 
magnetostatic SPS solution ([28]) as well as the require- 
ment of 3t{K ) > 1 constrain the parameter regime of 
this quasi-magnetostatic solution. In reference to equa- 
tions (6) and (7) of Lou & Wang (2007) corresponding to 
Q = here, we have explicitly discussed the parameter 
regime where two real roots K > 1 exist. We therefore 
expect that at least for a sufficiently small Q ^ 0, it 
would be still possible to have two real roots of K > 1 
given by quadratic equation ([33]) . Of course, it is also 
possible to have one root K > 1 and the other root 
K < 1. For example, for 7 = 1.4, n = 1.4 and h = 0, we 
have a positive K — 1.056 and a negative K = —1.199; 
with other parameters the same but h = 1, the two if 
roots are if = 1.0373 and —1.180, respectively. 



7 (3n - 2) 



1 1/(7-1) 



(7-1) 



m(0) 
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XX 



-(27-3)/(7-l) 



(37) 



Here m(0) is the value of the total reduced enclosed 
mass m(x) as x approaches 0, representing an increasing 
point mass at the very centre. This asymptotic solution 
at small x is valid for 3/2 < 7 < 5/3. The corresponding 
radial profile of the random magnetic field can be readily 
inferred from equation (|11| . 

3.1.5 Solutions of MHD Thermal Expansion at Large x 



The asymptotic MHD thermal fall solutions (136|) and 
(l37|) at small x in the preceding subsection 3.1.4 can be 
integrated numerically and connected with the asymp- 
totic MHD thermal expansion solution at large x (see 
Figure [6]). For this latter solution at large x, the reduced 
radial flow velocity v(x) becomes linear in x while the 
mass density approaches zero for x — > 00. The force that 
is dominant here is the thermal pressure gradient, driv- 
ing the gas into expansion. In this case of v ~ ex with 
c being a constant coefficient, we obtain from coupled 
nonlinear ODEs ([9]) and (fi"4|) at large x 



(38) 



(n — c)xa' — ca = —2(1 — c)a 



3.1.4 MHD Thermal Fall Solutions at Small x 

Lou & Li (2007 in preparation) report a possible asymp- 
totic thermal fall solution for small x without magnetic 
field; for this solution, the thermal pressure gradient force 
almost balances the gravitational force, yet the mass den- 
sity and radial infall speed still approach infinity at small 
X. In the case of asymptotic MHD thermal fall solutions, 
magnetic force is much smaller than both the pressure 
gradient force and the gravity at small x. The leading 
nonlinear ODEs for small x read 



vv' + P'/a - va/(3n - 2) = , 



(va)' + 2va/x 2 = 



(35) 



where the effect of magnetic field does not appear explic- 
itly (i.e., the absence of h parameter) in the regime of 
small x. This pair of coupled nonlinear ODEs (f34|) and 
(j35|) at small x are the same as those in Lou & Li (2007 
in preparation), and the corresponding asymptotic solu- 
tion is therefore the same to the leading order, namely 



7 (3n - 2) 
(7-1) 



-1/(7-1) 



m(0) 



-(9-1)7(7-1)^,-1/(7-1) 
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(n — l)ca; — (n — c)cx + /3 /a = 



(39) 



By further assuming the power-law form of a ~ Ex p at 
large x with E and P being two parameters and using 
equations ([8]) and (p~2|) . we obtain from equations ([38]) 
and 



(40) 



P = -(3g-2)/(l-n + 3ng/2) , 



E L 



VZnq / 2 (n-c) q {2 + P) =c(l-c) 



(34) and 



P = (3c-2)/(n-e) 



(41) 



(42) 



Parameters P, c and E are then readily determined by 
three algebraic equations PD|) — (j-4"2")) once the values of 
q and n parameters are specified. One can easily fix P 
from equation (|40|) first, then calculate c from equation 
([42]) and finally obtain E by equation ([41]) . This solution 
for MHD thermal expansion is valid when 7 > 4/3 (i.e., 
q > 2/3), because we require P < for a converging 
a(x) at large x. Again, the corresponding radial profile 
of magnetic field can be readily inferred from equation 

dnj. 

Asymptotic solution here should be compared with 
asymptotic solutions ([26]) and ([27]) at large x when n < 1 . 



3.2 MHD Solutions with Weak Thermal Pressure 

For MHD solutions involving weak thermal pressure, pa- 
rameter q should not enter the leading order terms of the 
asymptotic solution. 

3.2.1 MHD Free-Fail Solutions at Small x 

For asymptotic central free-fall solutions at small x, first 
found by Shu (1977) under the isothermal approxima- 
tion, the gravity force is virtually the only force in ac- 
tion, and the radial velocity and mass density profiles 
both diverge in the limit of x — » + . The leading order 
ODEs at small x then read 

2a a 2 



a = 



x (3n — 2)v 



(43) 



(3ti-2) ' 
and the corresponding asymptotic solution appears as 
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a{x) 



v(x) 



(3n - 2)m(0) 
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with an integration constant m(0), representing an in- 
creasing point mass at the very centre. The special isother- 
mal case of n — 1 and 7 = 1 was first studied by Shu 
(1977) in the context of star formation (Shu et al. 1987). 
It is particularly interesting to note that it is now possi- 
ble to have 7 = 1 and n > 2/3 case, still corresponding 
to a non-isothermal gas flow (i.e., q ^ and a variable 
sound speed). On the other hand, it is also possible to 
have n = 1 and 7 > 1 case with q ^ 0. With various 
combinations of sensible parameters, we can now readily 
construct such kind of solutions numerically. In terms 
of modeling the MHD processes of star or core forma- 
tion in magnetized clouds (Zhou et al. 1992; Shen & Lou 
2004; Fatuzzo et al. 2004; Lou & Gao 2006), our gen- 
eral polytropic MHD model framework is more versatile 
including the possible role of a completely random trans- 
verse magnetic field. In particular, we can model various 
molecular spectral line profiles in star forming regions; 
we shall pursue this exploration in separate papers. 

To the leading order, this asymptotic MHD free-fall 
solution does not involve polytropic index 7 and the cor- 
responding profile of magnetic field can be readily de- 
termined by equation (jTTJ) . In this case, magnetic field 
does not play a dynamically important role but may re- 
veal diagnostic features if shock accelerated relativistic 
electrons are present. We now consider the parameter 
regime that allows for this asymptotic MHD free-fall so- 
lution at small x. Substituting equations (|4"4"|) and (|4"S"]) 
into coupled nonlinear MHD ODEs (fl~5"|) and (fi"6l) and 



requiring the emergence of equation (|43j) for small x val- 
ues with consistent orders of magnitudes for the higher 
order terms, we obtain two inequalities 



n > 2/3 



and 



7< 5/3 



(46) 



these requirements appear to be the same as those of 
Suto & Silk (1988), but the polytropic equation of state 
adopted is different between theirs and ours, that is, it 
is no longer necessary to impose the condition n + 7 = 2 
or q = in our more general MHD polytropic model 
formulation. 

3.2.2 Strong- Field Asymptotic MHD Solutions in the 
Small-x Regime 

The strong-field asymptotic MHD solutions in the small- 
x regime (see Wang & Lou 2007 for the q = counter- 
part of this asymptotic solution in a conventional poly- 
tropic gas) are central or collapse solutions with the mag- 
netic Lorentz force against the self-gravity and with their 
strengths being comparable in magnitudes; both forces 
overwhelm the thermal pressure gradient force. For such 
strong-field asymptotic solutions, the velocity profile re- 
mains finite but the mass density profile becomes di- 
vergent as x — ► + . More specifically in the regime of 
x — > + , we obtain leading-order terms for the three co- 
efficients in equation (|16|) . namely 



X(x, a,v) ~ hax 2 



A(x, a,v) 



(nx-v) 2 



(3n - 2) 



a — 2ha x 



V(x, a, v) ~ 2hax 2 (l - n) + (na: ^ a , (47) 

(3n — 2) 

such that the two coupled nonlinear MHD ODEs in (fT5|) 
can be simplified to 



v(x) 



and 



(3n - 2) 



h±{h 2 - Ah) 



1/2 



a(x) = DoX -5/2TV^IK/{2h) 



(48) 



(49) 



to the leading order of small x, where Dq is an integration 
constant; the corresponding magnetic field strength pro- 
file is simply characterized by equation (jlip . For strong- 
field asymptotic MHD solutions (05]) and (@3]) to be phys- 
ically valid, it is necessary to require that h > 4, corre- 
sponding to a sufficiently strong magnetic field regime 
and hence the name of this type of asymptotic MHD 
solutions. In a wide range of astrophysical systems, the 
typical value of h ranges from ~ 10~ 3 to ~ 10 5 . The 
strong-field situations may happen for magnetosphcric 
plasmas surrounding magnetic white dwarfs, radio pul- 
sars, anomalous X-ray pulsars (AXPs), magnetars and so 
forth. Typical magnetic field strengths are ~ 10 7 — 10 9 G 
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for magnetic white dwarfs, ~ 10 9 — 10 10 G for millisec- 
ond radio pulsars, ~ 10 11 — 10 12 G for radio pulsars, 
- 10 12 - 10 13 G for AXPs, and - 10 13 - 10 15 G for mag- 
netars. 

To the leading order, the polytropic index 7 does not 
appear in the strong-field asymptotic solution, because 
the thermal pressure force is much weaker than both 
gravity and magnetic force. By requiring [3 — rrfia 1 <C 
a 2 x 2 in the derivation, the parameter regime for these 
strong-field asymptotic MHD solutions at small x to be 
valid is 




5 (h 2 - Ah) 1 / 2 
2 T 2h 



3nq 

~2 



-n+— -)+3q-2>0 . (50) 



Fig. 1 Semi-complete MHD solutions with inner free-fall 
asymptotic solutions and without MHD shocks. All these so- 
lutions involve a completely random magnetic field and a gen- 
eral polytropic gas. The relevant parameters are 7 — 1.1, n = 
0.85, q = —0.182, and h = 0.3. The two perpendicular dot- 
ted straight lines are abscissa and ordinate axes, respectively. 
The dash-dotted curves are the magnetosonic critical curves. 
In the upper panel, the MHD solution labelled with 'Shul' is 
a polytropic MHD counterpart of the isothermal expansion- 
wave collapse solution (EWCS; Shu 1977), and is co nstru cted 
by integrating from large x asymptotic solutions (I26|l and 
(|27[) with parameters A slightly larger than Ao = 1.842 and 
B = 0; the MHD solution labelled with 'Shu2' is constructed 
in a si mil ar m anne r for A = 2 and B = in asymptotic solu- 
tions (J26J and l[2T|); the MHD solution labelled with 'Hunterl' 
is constructed with two parameters A — 2 and B — —0.3; 
the 'Hunter2' solution is constructed with A = 1.6 < Aq 
and B = —1; and the 'Hunter3' solution with A = 3 and 
B = 1. We refer to these MHD counterparts as Hunter 
type solutions, because Hunter (1977) implicitly constructed 
such isothermal solutions without magnetic field in the com- 
plete space. In the lower panel, the MHD solution labelled 
with 'LSI' is constructed by integrating from (x,v,a) = 
(0.03092, —1.581, 56.94) on the magnetosonic critical curve 
inwards to small x and outwards to xf = 0.3, and by integrat- 
ing from (as, V, a) = (2.683, 0.9359, 0.4795) on the magne- 
tosonic critical curve inwards to the same if and outwards to 
large x; this solution goes smoothly across the magnetosonic 
critical curve twice. The other MHD solution labelled with 
'LS2' is constructed by integrating from (x, v, a) — (3.631 x 
10~ 6 , —4.944, 1.497x 10 6 ) on the magnetosonic critical curve 
inwards to small x and outwards to xf = 0.3, and by integrat- 
ing from (x,v, a) = (0.4357, -0.8466, 2.8839) on the magne- 
tosonic critical curve inwards to the same xf and outwards 
to large x (Lou & Shen 2004). The values of m(0) for the 
'Shul', 'Shu2', 'Hunterl','Hunter2','Hunter3', 'LSI' and 'LS2' 
MHD solutions are 1.286,1.565,1.777,1.767,1.864,0.08699 
and 9.763 x 10" 5 , respectively. The 'LSI' MHD solution has 
its large x asymptotic solution with parameter s A = 4.241 
and B = 1.843 in asymptotic solution (|26[) and (|27|) ; and the 
'LS2' MHD solution has its large x asymptotic solut ion with 
A = 0.4973 and B - -2.431 in asymptotic solution ([26) and 
((27)) . As examples of illustration, the 'LSI' and 'LS2' MHD so- 
lutions cross the magnetosonic critical curve twice smoothly. 



By setting q = and thus n + 7 = 2, inequality (fSTJ)) 
reduces to inequality (69) of Wang & Lou (2007) for a 
conventional polytropic gas permeated with a completely 
random transverse magnetic field. 

Numerical examples of using strong-field asymptotic 
solutions (|48p and (|4"9")l with the upper signs are shown in 
Figure^ These solutions (labelled by 'WL1' and 'WL2') 
pass through the magnetosonic curve smoothly and join 
the asymptotic solutions (|26|) and (|27|) at large x. One 
can also construct MHD shock solutions in this scheme 
numerically. 

We note empirically that numerical integrations from 
very small x outwards using strong-field asymptotic solu- 
tions (|48|) and (|49)l have a tendency to be unstable, espe- 
cially for the lower-sign solution. On the other hand, we 
are unable to obtain the lower-sign solutions in (|48[) and 
(l4"9"|) by numerically integrating from the magnetosonic 
critical curve to small x so far. This problem might be re- 
lated to properties of perturbations in a self-similar flow 
(Lou & Bai 2008 in preparation). 

Lou & Wang (2007) proposed a class of self-similar 
MHD rebound shock models for supernovae based on 
quasi-static MHD shock solutions with q — and fo- 
cussed on the origin(s) of intense magnetic fields on rem- 
nant compact objects such as radio pulsars (magnetic 
field strengths of ~ 10 11 — 10 12 G) and magnetic white 
dwarfs (surface magnetic field strengths of ~ 10 7 — 10 9 G). 
It can be possible that MHD rebound shock models based 
on our strong-field solutions (see Figure [5]) with g ^ 
are relevant to strong magnetic fields of ~ 10 13 — 10 15 G 
observationally estimated for magnetars and anomalous 
X-ray pulsars (AXPs). This would correspond to strong 
surface magnetic field strengths higher than several thou- 
sand gauss on progenitor stars (e.g., magnetic Ap stars). 
We shall discuss this interesting problem in separate pa- 
pers. 



4 Global Semi-Complete MHD Solutions 
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and Astrophysical Applications 

4.1 MHD Expansion- Wave Collapse Solutions 

Shu (1977) constructed the expansion-wave collapse so- 
lution (EWCS) for a self-similar isothermal gas flow and 
later developed the so-called inside-out collapse scenario 
for protostar formation (see Shu, Adams & Lizano 1987 
and extensive references therein). Cheng (1978) presented 
the polytropic generalization of EWCS. Chiueh & Chou 
(1994) studied the MHD generalization of an isothermal 
EWCS. In reference to the work of Yu & Lou (2005) 
in this context, we can now construct the general poly- 
tropic MHD EWCS counterpart here and show a specific 
example labelled with 'Shul' in Figure [TJ In our general 
polytropic MHD EWCS, the outer portion is the outer 
part of a magnetized SPS (no singularity is actually in- 
volved because x — is excluded), while the inner por- 
tion approaches an MHD asymptotic free-fall solution 
with the random magnetic field being advected radially 
inward. From the perspective of theoretical model de- 
velopment, we can now test thermodynamic properties 
of magnetized clouds by comparing with observational 
inferences. Using the isothermal version of this solution, 
Shu (1977) outlined a physical scenario for protostar for- 
mation, which can now be further extended on the ba- 
sis of our generalized polytropic MHD model solution. 
At the beginning, the magnetized SPS is static; then at 
t = 0, a disturbance may take place and a magnetized 
core collapse under self-gravity occurs in an inside-out 
manner; the 'bounding surface' (or the stagnation sur- 
face) of the magnetized collapsing sphere travels outward 
at the magnetosonic speed (not a constant in general) in 
a self-similar fashion. Shu (1977) utilized the isothermal 
version of this EWCS solution to model the core collapse 
of a molecular cloud in protostar formation. Zhou et 
al. (1993) attempted to test Shu's isothermal dynamical 
model by fitting the observed molecular spectral line pro- 
files. However, it should be noted that by incorporating 
a non-isothermal temperature profile inferred from data, 
this isothermal model test by fitting molecular spectral 
line profiles is not fully self-consistent. In contrast, with a 
general polytropic model and the ideal gas law, the tem- 
perature variation is a natural consequence. It would be 
highly desirable to carry out a parallel analysis similar to 
that of Zhou et al. (1993) but for a general polytropic gas 
cloud model so that profiles of temperature, density and 
flow speed are all self-consistently prescribed. By setting 
h = for the absence of magnetic field, our general- 
ized polytropic EWCSs are the same as those of Cheng 
(1978). The outer magnetized SPS involves a completely 
random magnetic field and thus a global quasi-spherical 
symmetry. As a result, small-scale random flows with a 
zero mean are expected. In this sense, the outer magne- 
tized SPS portion here provides a large-scale mean pro- 
file. The presence of a magnetic field would make a cloud 
to behave more fluid like on large scales (the solar wind 



is an example) and can give rise to radiative diagnostics 
when relativistic electrons are involved. 

The main difference between our general polytropic 
MHD EWCS here and former EWCSs lies in the degree 
of freedom. For every determined MHD flow with a set 
of parameters {7, n, h}, there is at most one EWCS and 
these three parameters are now allowed to change inde- 
pendently. In the isothermal EWCS case of Shu (1977), 
we have 7 = 1, n = 1 and ft, = 0. In the conventional 
polytropic EWCS case of Lou & Wang (2006), we have 
7 + 11 = 2 (q = 0) and h = 0. In the conventional poly- 
tropic MHD problem studied by Wang & Lou (2007), 
q = or n + 7 = 2 is required and there are two param- 
eters 7 and h that are allowed to vary independently. 
In the present generalized polytropic MHD model, the 
parameter n is also allowed to change and is the key ex- 
tension of this paper for this model problem. Therefore, 
all polytropic MHD solutions in this analysis have one 
more degree of freedom for fitting observational data. 

4.2 Inner Free-Fail with Outer Inflow/Outflow 

Shu (1977) also constructed global semi-complete isother- 
mal solutions that approach a central free fall at small 
x, but with B = 0, h = 0, n = 7 = 1 (or q = 0) 
in asymptotic solution (|27|) and here A > Ao for the 
outer asymptotic singular isothermal sphere (SIS) solu- 
tion (n.b., the singularity at x — is actually excluded 
in this construction). In fact, Shu used this sequence of 
solutions to introduce the EWCS as the limiting case 
and to suggest an inside-out collapse scenario for star 
formation. As an example, we show a generalized poly- 
tropic MHD extension of such solutions in Figure [IJ la- 
belled by 'Shu2' with B = and A > A ; it is possible 
and straightforward to construct a sequence of such so- 
lutions with B = by choosing different A values larger 
than Ao. In addition, there is still the case of B ^ in 
asymptotic solution l|27[). which is qualitatively similar 
yet has a constant flow speed at large x in the isothermal 
case for outer envelopes (e.g., Lou & Shen 2004; Shen & 
Lou 2004; Fatuzzo et al. 2004; Yu & Lou 2005; Yu et al. 
2006). There are increasing observational evidence indi- 
cating that star forming clouds do have systematic flows 
far away from the core (e.g., Fatuzzo et al. 2004) and our 
model can accommodate either inflows or outflows. 

In general polytropic cases with n < 1, the radial 
flow speed remains finite and approaches zero at large x. 
With n = 1, the radial flow speed remains finite and ap- 
proaches a constant value B at large x. While for n > 1, 
the asymptotic flow term associated with this B coeffi- 
cient in asymptotic solution (|27p diverges at large x; to 
model a real system, one then needs to set B equal to 
zero. This parameter B represents a component of ra- 
dial flow speed which satisfies the momentum equation 
as if there is no force in action. Because Hunter (1977) 
first constructed global semi-complete isothermal solu- 
tions with this behaviour, we label 'Hunterl', 'Hunter2' 
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Fig. 2 Semi-complete MHD shock solutions with free falls as 
x — > with four parameters 7 = 1.05, n = 0.9, q = —0.143, 
and h = 0.1. In the upper panel, the MHD shock solution 
labelled with 'BL1' is constructed by integrating from a 
magnetosonic critical point (x, v, a) — (0.5, —0.6520, 3.053) 
inwards to small x and outwards to a MHD shock point 
(x S 2,v S 2,a S 2) — (0.7, — 0.2 731, 2.19 1) downstream; using 
MHD shock conditions (|22[) — (|24[) . we obtain the corre- 
sponding upstream physical variables and thus integrate 
from (x sl ,v s i,a s i) = (0.70015,-0.6761,1.515) outwards to 
large x. The 'BL2' MHD shock solution is constructed by 
the same inner portion with a different MHD shock point at 
(x S 2,v S 2,a S 2) = (1.9,0.8275,0.8481); again using MHD shock 
conditions (|22[) — (|24[) . we determine the upstream physical 
parameters at (x s i,v s i,a s i) = (1.90061,0.4034,0.5728) 
and then integrate outwards to large x. The 'BL1' MHD 
solution has a large x asymptotic sol ution with A = 0.6536 
and B = —1.780 in expressions ((26} and (|27[1 ; while these 
parameters for 'BL2' MHD sho ck s olutio n a re A — 2.168 
and B = 0.4534 in expressions (|26p and (|27p . In the lower 
panel, the 'TH1' and 'TH2' MHD solutions both have the 
same outer SPS portion with v — and a = Aox~ 2 ' n 
as described by equation 1)280 ; this is the outer part of a 
magnetized SPS with t he r adial profile of a random magnetic 
field given by relation (jXTJ . The MHD solution labelled with 
'TH1' has upstream and downstream MHD shock points 
(x sl ,v sl ,a sl ) = (1.7555,0,0.5311) and (x S 2,v S 2, a s2 ) = 
(1.7514,0.8358, 1.130), and crosses the magnetosonic critical 
curve at (x,v,a) = (0.1376,-1.081,12.20). The MHD 
solution labelled with 'TH2' has upstream and downstream 
MHD shock points (x s i,v s i,a s i) = (1.4766,0,0.7800) 
and {x s2 ,v S 2,a S 2) = (1.4760,0.4623,0.1963), and 
crosses the magnetosonic critical curve at (x,v,a) = 
(7.825 x 10" 6 , -1.081, 4.972 x 10 5 ). The 'BL3' solution has 
shock points (x sl ,v sl ,a sl ) = (0.4758,-1.298,3.289) 
and (x 3 2,v 3 2,a 3 2) = (0.4743,-0.2907,7.889), and 
crosses the magnetosonic critical curve at (x,v,a) = 
(4.885 x 10" 4 , -1.945, 4.474 x 10 3 ). The values of m(0) for 
solutions 'BL1', 'BL2', TH1', 'TH2' and 'BL3' are 0.6987, 
0.6987, 0.2665, 8.461 xl0~ 5 and 2.589 x 10" 3 . The values of 
shock parameter A for these solutions are 1.000210, 1.000320, 
1.00232, 1.000390 and 1.00303, respectively. 



and 'Hunter3' in Fig. Q]as examples of this type of flow 
solutions. Here the 'Hunter2' solution is one with a less 
dense envelope than the SPS yet with an inflow speed; 
the 'Hunter3' solution is one with a more dense envelope 
yet with an outflow speed. Both involve free-fall solutions 
in the central core collapse region. 

Qualitatively, this type of general polytropic MHD 
solution has either an inflow or an outflow in the outer 
envelope and a free-fall inner core (Lou & Shen 2004; 
Shcn & Lou 2004; Yu & Lou 2005; Yu et al. 2006; Lou 
& Gao 2006). As the presence of inflows is inferred ob- 
servationally in star-forming regions and molecular cloud 
cores (see, e.g., Fatuzzo, Adams & Myers 2004 and refer- 
ences therein), our model analysis here is more general in 
the following four aspects, namely, (i) the inclusion of a 
completely random magnetic field, (ii) the construction 
of an MHD shock, (hi) the possibility of a constant 'spe- 
cific entropy' everywhere at all time, and (iv) the possi- 
bility of 'specific entropy' conservation along streamlines 
yet with a variable 'specific entropy' in space and time. 
Physically, unless for a 'black hole' at the centre, the 
central free-fall of gas will ultimately lead to a 'sphere of 
transient activities' around the core region; the change of 
equation of state and/or the ignition of thermal nuclear 
reaction mark the onset of a star formation process. In 
other words, we expect that the mathematical singular- 
ity of a free-fall MHD solution as x — ► + will be taken 
care of by other relevant physical processes when a proto 
star forms. 

Apart from the index parameters n and 7 that are 
allowed to vary independently for different polytropic 
MHD flows, there are still two parameters A and B in 
asymptotic MHD solutions (|2"6")) and (|2~T|) . In a contin- 
uous manner, these two parameters can change and be 
mapped to a central mass point m(0) or central mass ac- 
cretion rate. This mapping can be readily determined by 
numerical integration. Fatuzzo et al. (2004) made a semi- 
analytical estimate of this mapping. In general, larger A 
values or smaller B values will lead to a larger m(0) as 
expected on intuitive ground. Observationally, this cor- 
responds to various central mass accretion rates in pro- 
tostar forming clouds. 

4.3 Envelope Expansion with Core 
Collapse in Magnetized Clouds 

Lou & Shen (2004) constructed isothermal EECC solu- 
tions using a solution matching technique in the a — v 
phase diagram (see also Bian & Lou 2005 and Wang & 
Lou 2007). We constructed such kind of general poly- 
tropic MHD solutions and present the first two of this 
series of MHD solutions in the lower panel of Figure [1] 
The first EECC solution labelled with 'LSI' is a MHD 
solution with an inner free fall (i.e., magnetized core col- 
lapse) and an outer outflow (i.e., magnetized envelope ex- 
pansion) . This kind of evolution behaviour may be qual- 
itatively applicable in asymptotic giant branch (AGB) 



stars, post-AGB stars, and/or protoplanetary nebulae 
(PPNe); along this line, Lou & Shen (2004) proposed 
to utilize EECC solutions (including the second 'LS2' 
solution which is qualitatively different) to grossly catch 
key features of MHD collapses and flows. Also in the 
context of a collapsing magnetized cloud to form a pro- 
tostellar core, general polytropic MHD EECC solutions 
with o^O exemplified here appear more general to ac- 
count for various plausible situations (Shen & Lou 2004; 
Lou & Gao 2006). In all these considerations, the inner 
singularity as x — > + can be removed once a departure 
from self-similarity occurs and/or a different equation of 
state is adopted. 

Except for the additional scaling parameter n, for 
every fixed MHD flow profile there exists at most a series 
of discrete solutions, each with a qualitatively different 
manner (e.g., different times to cross the v = axis). It 
is common to this kind of flows (see also Hunter 1977; 
Lou & Shen 2004; Bian & Lou 2005) that the a versus 
v phase diagram obtained by integrating from the inner 
portion (from x values smaller than the meeting point 
xf) tends to spiral into a specific point. This point is on 
the phase curve obtained by integrating from the outer 
portion (from x values larger than the meeting point Xf)- 



4.4 Expansion of MHD Shocks into an Outer Static 
Envelope Configuration 

Tsai & Hsu (1995) constructed an isothermal self-similar 
shock solution which connects an outer static configura- 
tion to a central free-fall solution with a shock. Bian & 
Lou (2005) further constructed various possible isother- 
mal shock solutions in a more comprehensive manner. 
We present here the first three of such MHD shock solu- 
tions in the lower panel of Figure [2] Based on extensive 
numerical exploration, Tsai & Hsu (1995) also suggested 
that in a star-forming cloud, when a strong burst of ther- 
mal nuclear energy is released from the core, instead of 
a smooth evolution in the form of EWCS as discussed 
by Shu (1977), a shock can gradually emerge and travels 
outward in a self-similar manner. Due to a sharp density 
profile of the singular isothermal sphere (the counterpart 
of which is SPS in our model analysis here), the shock 
travels long and evolves into a self-similar shock. They 
provided numerical simulations to support this scenario 
of shock initiation and evolution. 

The degree of freedom of this series of generalized 
polytropic MHD shock solutions is like that of the MHD 
EECC solutions. Because the outer SPS is prescribed, 
the MHD shock solutions need to be properly matched in 
the a versus v phase diagram. Similar spiral-in features 
in the a versus v phase diagram (e.g., Lou & Shen 2004) 
may exist, leading to a series of discrete MHD shock 
solutions. Corresponding to each of such MHD shock so- 
lutions, there is a specific central mass accretion rate. 
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Fig. 3 Semi-complete self-similar polytropic solutions 
with twin MHD shocks for parameters 7 = 1.05, n = 0.9, 
q — —0.143, h — 0.1. The dash-dotted curve is the 
magnetosonic critical curve. The inner portion is the 
same for the two MHD solutions labelled with 'BL4' and 
'BL5', i.e., the MHD shock points for the inner shock are 
(x S 2,v s2 ,a S 2) = (0.4129,-0.4458,6.949) for the downstream 
side and (x sl ,v sl ,a sl ) = (0.4133,-1.147,3.743) for the 
upstream side, and both MHD solutions also smoothly cross 
the magnetosonic critical curve twice at the same two points 
at (x,v,a) = (1.425 x 10" 4 , -2.164, 2.126 x 10 4 ) and at 
(x,v,a) = (0.7,-0.4543,2.146). The common inner portion 
of both MHD solutions has two nodes (i.e., v = 0) for x < 0.3 
and diverges as 1 —> 0. Parameter A for the inner MHD 
shock is 1.00100 and parameter m(0) for the central free-fall 
solution is 9.338 x 10 -4 . The MHD solution labelled with 
'BL4' has the outer MHD shock at (x S 2, v S 2,a 3 2) = 
(1,-0.03170,1.512) for the downstream side and 
(x s i,v sl ,a Bl ) = (1.0000961, -0.3316, 1.144) for the upstream 
side; while the 'BL5' MHD solution has the outer MHD shock 
at (x S 2,v 3 2, ct a 2) — (1.9,0.7693,0.8311) for the downstream 
side and (x sl ,v sl ,a sl ) = (1.90019,0.4803,0.6358) for the 
upstream side. The values of A for the two outer MHD 
shocks of 'BL4' and 'BL5' are 1.0000961 and 1.000102. The 
'BL4' MHD solution has a large x asymptotic solutions (|26[) 
and lO with A = 1.076 and B = -0.9545 (inflow), while for 
the 'BL5' MHD solution, the two corresponding parameters 
are A = 2.403 and B = 0.6354 (outflow). 



4.5 Generalized Polytropic MHD Shock 
Solutions with Inner Free Fall 

Bian & Lou (2005) explored various isothermal hydro- 
dynamic shock solutions and noted astrophysical appli- 
cations to AGB, PNe, protostar formation, quasars and 
supernova explosions and so forth. We here present two 
generalized polytropic MHD shock solutions 'BL1' and 
'BL2' as illustrative examples in Figure [U Within the 
framework of our model analysis, it is possible to con- 
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struct a variety of general polytropic MHD shocks adapted 
to various astrophysical flow situations. 

On the basis of a conventional polytropic hydrody- 
namic shock model with n+"f = 2 and thus q — 0, Lou & 
Gao (2006) examined observationally inferred informa- 
tion for star-forming cloud cores. They noted that star- 
forming regions are well studied in the inner core and 
outer edge regions, which may be utilized to constrain 
or test self-similar core collapse scenario with or without 
inflows and/or outflows (Shen & Lou 2004; Fatuzzo et al. 
2004). While globally smooth solutions such as EECC so- 
lutions have less degrees of freedom when they encounter 
the sonic critical curve (which occurs often), the possi- 
ble existence of shocks is physically plausible and math- 
ematically convenient to join different inner and outer 
asymptotic solutions into a global shock flow solution 
(i.e., by crossing the sonic critical curve with a shock). 
Lou & Gao (2006) also compared polytropic model with 
observations. By taking into account of radiative trans- 
fer, we shall further fit various observed molecular spec- 
tral line profiles for given underlying general polytropic 
MHD shock flows with q ^ 0. Thus shock solutions pro- 
vide a simple and direct model scenario for star-forming 
cloud cores that appear grossly quasi-spherical on large 
scales. 

Once parameters 7, n and h are prescribed, a global 
semi-complete MHD shock solution has two degrees of 
freedom, namely, one can choose the location where the 
solution smoothly crosses the magnetosonic critical curve, 
and one can also choose the MHD shock point; this is il- 
lustrated in the two panels of Figure [2] where two MHD 
shock solutions have different outer envelopes. The loca- 
tion of smoothly crossing the magnetosonic critical curve 
determines the parameter m(0) (related to the central 
mass point and mass accretion rate) in the core and the 
complete choice corresponds to a specific set of param- 
eters A and B for the outer large x asymptotic MHD 
solution (J26J) and (|2?)), The plethora of this type of MHD 
shock solutions can accommodate various astrophysical 
situations including protostar-forming cloud cores. 



4.6 Twin Polytropic MHD Shock Solutions 

Bian & Lou (2005) also constructed isothermal hydro- 
dynamic twin shock models, in which two shocks ap- 
pear in a self-similar flow and the solutions cross the 
sonic critical curve thrice, namely, once smoothly and 
twice with shocks. For general polytropic MHD flows, we 
present the counterpart solutions in Figure OS This type 
of general polytropic MHD shock solutions further ex- 
pand the solution space, providing more plausible mod- 
els of quasi-spherical processes involving random mag- 
netic field, gravitational core collapse and far-away in- 
flows/outflows. We note that this general polytropic MHD 
twin shock model differs from the so-called forward-reverse 
shock pair (see, e.g., Chevalier 1982), because both shocks 



here are 'forward shocks' in the sense that the shock 
moves forward relative to the local MHD flow. Concep- 
tually, we emphasize the possibility of twin or multiple 
general polytropic MHD shocks in a magnetized flow sys- 
tem with q y^ (see Yu et al. 2006 for the isothermal 
case). 

When three parameters 7, n and h are chosen, one 
can still adjust the place where the MHD flow crosses the 
magnetosonic critical curve smoothly, and then choose 
the place where an outer MHD shock appears; the for- 
mer choice gives the central mass point m(0) or central 
mass accretion rate, while for the latter choice, the shock 
location x also corresponds to the outward travel speed 
of this MHD shock. Figure [3] shows that different choices 
of outer MHD shocks correspond to very different outer 
MHD flow profiles in the envelope. 



4.7 Inner Quasi-Magnetostatic 
Solutions with Outer Inflows 

Using a conventional polytropic hydrodynamic formula- 
tion with q — 0, Lou & Wang (2006) first presented the 
quasi-static asymptotic solution in the regime of small 
x and then constructed a rebound shock model to catch 
certain gross features of a class of supernovae. We at- 
tempted to relate the mass of a progenitor star and 
the type of a remnant compact object (such as a white 
dwarf, a neutron star or a black hole) left behind af- 
ter a gravitational core collapse and a subsequent emer- 
gence of rebound shock. The above model can be ex- 
tended to include a completely random transverse mag- 
netic field with quasi-spherical symmetry (Lou & Wang 
2007). In this MHD rebound shock model, we explore 
the origin(s) of intense magnetic fields on compact ob- 
jects from the perspective of fossil field associated with 
massive progenitor stars. As an important part of fur- 
ther model development, we have also constructed gen- 
eralized polytropic MHD solutions that smoothly cross 
the magnetosonic critical curve and approach this quasi- 
magnetostatic asymptotic solution (|29|) as shown in Fig- 
ure [4] and labelled with 'LWF. In the regime of small x 
and with complex coefficients K and L, the MHD flow 
velocity actually oscillates with a decreasing magnitude 
as x — > + (not easily seen here; see Lou & Wang 2006 
for a similar example in detail). Qualitatively speaking, 
this generaized polytropic MHD solution represents a 
far-away or initial inflow leading to an eventual quasi- 
magnetostatic inner core and carries a desirable feature 
of forming a central magnetostatic configuration instead 
of a divergent free fall towards a central mass point. 
Therefore, this type of generalized polytropic MHD so- 
lutions should also provide a physically plausible model 
framework to study protostar formation processes (Shen 
& Lou 2004; Lou & Gao 2006); we shall discuss this im- 
portant application in separate papers. 

This type of general polytropic MHD solution has 
one degree of freedom in construction, namely, the choice 



of the point where a solution crosses the magnetosoni 
critical curve. 



4.8 Inner Quasi-Magnetostatic Solutions 
with Rebound MHD Shocks 

Lou & Wang (2006) and subsequently Lou & Wang (200' 
proposed a conventional polytropic gas model (i.e., q = 
and thus n + 7 = 2) for rebound shock process in 
class of supernova explosions. In particular, Lou & Wan 
(2007) have incorporated a random magnetic field (e.g 
Zel'dovich & Novikov 1971) and attempted to addres 
the fundamental issue on the physical origin of stron 
magnetic field of a compact object left behind a reboun 
MHD shock. Based on our estimates, we proposed tha 
for a progenitor star of initial mass in the range of <■ 
6 — 8M0, an MHD rebound shock initiated during th 
core collapse may eventually leave behind a magneti 
white dwarf with magnetic field strengths in the rang 
of ~ 10 6 — 10 8 G depending on surface magnetic fieL 
strengths of the progenitor star. Similarly, for a pre 
genitor star of initial mass in the range of > 8M Q , a 
stronger MHD rebound shock initiated during the grav- 
itational core collapse may ultimately leave behind a 
neutron star with magnetic field strengths in the range 
of ~ 10 11 — 10 12 G depending on surface magnetic field 
strengths of progenitor stars. In reference to the key re- 
sult of Lou & Wang (2007), the generalized polytropic 
MHD model here allows n + 7 7^ 2 and thus g ^ 0. 
The major physical consequence is that 'specific entropy' 
is conserved along streamlines but 'specific entropy' is 
not uniformly distributed in space and time during the 
gravity-induced core collapse and rebounce within a pro- 
genitor star. As an example of illustration, we display 
our generalized polytropic rebound MHD shock solution 
labelled with 'LW2' in Figure 5) This solution carries 
a feature that a strong MHD shock emerges surround- 
ing an inner magnetized core, ploughs through an in- 
falling magnetized outer envelope (either with a stel- 
lar wind or just with an inflow) and leaves behind a 
quasi-magnetostatic compact configuration of high den- 
sity. The 'rebound' MHD shock here is a physical sim- 
plification of a neutrino-driven shock, as opposed to the 
earlier concept of a 'prompt' shock. One important flexi- 
bility here is that the constraint n + 7 = 2 is now relaxed 
such that relevant coefficients allow for various plausible 
combinations. In the analysis of Lou & Wang (2007), 
we have already pointed out that the asymptotic scal- 
ing laws of the mass density profile and magnetic field 
profile in terms of radius r is independent of the equa- 
tion of state, which is apparent from the analysis of this 
paper. This is the crucial difference between this model 
under the present formulation and those of former anal- 
yses. For example, in our former model (Lou & Wang 
2007), it appears that 
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Fig. 4 Semi-complet e so lutions with quasi-magnetostatic 
asymptotic solutions (|29p at small x for 7 = 1.2, n = 
0.7, q = -2 and h = 1. We then have Q = -3.308 
by definition (|32p and Aq = 1.078 in expression (|29|l . 
The dash-dotted curve is the magnetosonic critical curve. 
In this c ase, the parameter K obtained from qu adra tic 
equation (1331) for the quasi-magnetostatic solution (|29p is 
K = Ki + iK 2 = 1.357 + 0.834K. The MHD solution la- 
belled with 'LW1' is constructed by starting from (x,V,a) = 
(1, —1.016, 0.3583) on the magnetosonic critical curve and in- 
tegrating both inwards and outwards numerically. For the 
inner quasi-magnetostatic asymptotic solution at small x, 
parameter L is L - -0.02262 - 0.8951i, and thus N - 
—0.2365 — 1.459i; while for the outer large x MHD asymp- 
totic solution (|26p and (|27p . we have corresponding param- 
eters A = 0.7662 and B = -3.314 for an inflow. The 
polytropic MHD shock solution labelled with 'LW2' de- 
cribes a rebound MHD shock in a self-similar evolution (see 
Lou & Wang 2007 for more details); this solution is con- 
structed as follows. We start by integrating from (x, v, a) = 
(3, 0.8699, 0.1134) on the magnetosonic critical curve in- 
wards and then stop at (x S 2, v S 2,a S 2) = (1.8, 0.6877, 0.5869) 
which is regarded as the downstream side of an MHD shock 
and this segment of sol utio n is itself ignored. Using MHD 
shock conditions (I22p ~ (1241) . we then derive physical vari- 
ables (Xsi,Vsi,ct s i) = (2.623,-3.275,0.09576) upstream of 
the shock. Finally, we integrate inwards from the down- 
stream side with (x S 2,v S 2,a S 2) — (1.8, 0.6877, 0.5869) 
and outwards from the upstream side with (x s i,v 3 i,a s i) = 
(2.623, -3.275,0.09576). The parameter A for this polytropic 
MHD rebound shock is A = 1.457 and the outer portion 
has a large x asymptotic MHD solution with parameters 
A — 3.890 and B — —1.525 for an inflow. The inner quasi- 
magnetostatic solution has parameter L — 0.9797 + 0.4042i 
and thus N - 1.692 + 0.4427i Both MHD solutions oscillate 
in the small x regime in a self-similar manner (Lou & Wang 
2006; Wang & Lou 2007). 
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and we need n to approach 2/3 for the strongest mag 
netic field amplification or the fastest field variation fo: 
the stellar surface to the central core; this appears to con 
strain the value of polytropic index 7 in our former anal 
ysis with q = 0. However, in the current model frame 
work with 5^0, variations of n and 7 are no longe 
constrained in this regard. We can still let n approache 
2/3 but the 7 value is not necessarily close to 4/3. 

In constructing such general polytropic rebound MH] 
shock solution, we can choose the point where the solu 
tion would cross the magnetosonic critical curve smoothl 
(see Fig. @] caption for the solution construction pro 
cedure; the solution portion between the magnetosonii 
critical point and the downstream MHD shock point ii 
actually ignored later) and the MHD shock point. There 
fore except for the choice of parameters 7, n and h, wi 
still have two degrees of freedom for constructing sucl 
type of semi-complete solutions with q =/= 0. 

To describe a collapsing stellar core prior to the emer 
gence of a rebound shock for a supernova explosion, Gol 
dreich & Weber (1980) derived an exact homologous so 
lution using a polytropic model with 7 = 4/3 and con 
stant specific entropy everywhere within the stellar core, 
compared with numerical simulations (Van Riper & Ar- 
nett 1978; Bethe et al. 1979), and identified parameter 
range allowing for a homologous core collapse. In the 
light of current model analysis and complementary to 
the analysis of Goldreich & Weber (1980), Lou & Cao 
(2008) realized that it is not necessary to impose the 
condition of a constant specific entropy throughout the 
stellar core (also unlikely in reality) and one can still de- 
rive homologous solutions for a collapsing stellar core. 
Besides, we further obtain a broad class of self-similar 
solutions for 7 = 4/3 with or without shocks. 

4.9 Inner Magnetoaccretion with Outer Inflows 

Wang & Lou (2007) report a novel magnetoaccretion so- 
lution for a conventional polytropic gas (q = 0) with a 
completely random magnetic field. The isothermal MHD 
counterpart of this magnetoaccretion solution was de- 
scribed in Yu et al. (2006). We here construct the coun- 
terpart solutions (|4"5f and (|4*9")1 with the upper signs in 
a generalized polytropic (q ^ 0) MHD formulation and 
display two examples of these results in Figure These 
solutions involve fairly strong magnetic field (i.e., h > 4). 
Especially at small x, we would expect a 'sphere' of 
magnetic activities and transients within which the self- 
similarity and quasi-spherical symmetry are destroyed. 
However, the physical scenario corresponding to this so- 
lution, i.e., the possibility that the gravitational energy 
can be effectively converted into magnetic energy via 
magnetoacretion processes, is tantalizing. We shall fur- 
ther analyze this solution and its astrophysical applica- 
tions in separate papers. 

In constructing this type of general polytropic MHD 
solution, we can choose the point where the solution 
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Fig. 5 Semi-complete polytropic MHD sol utions wit h inner 
magnetoaccretion asymptotic solution (|48[) and (|49[) in the 
strong-field regime with 7 = 1.2, n = 0.7, q — —2 and 
h — 10. The dash-dotted curve is the magnetosonic criti- 
cal curve. For these parameters, the value of the prop ortional 
coefficient of v versus x is —0.1873 in solution (148 \ for the 
upper plus sign and this value is adopted for both WL1 and 
WL2 solutions. The MHD solution labelled with 'WL1' is con- 
structed by integrating from (x,v,a) = (1,-1.574,0.2671) 
on the magnetosonic critical curv e b oth inwards and out- 
wards. Parameter Do in solution (|49[1 is about 1.4 and the 
outer portion has large x asymptotic solution (|26[) and (|27[) 
with parameters A = 0.7301 and B = —5.155 (inflow). The 
'WL2' solution is constructed similarly, but the starting point 
is (x,v,a) = (1.5,-1.151,0.1452) o n th e magnetosonic criti- 
cal curve; parameter Do in solution (149 II is ab out 1.5, and the 
large x MHD asymptotic solution (|26p and (|27|) has param- 
eters A - 0.8738 and B = -3.494 (inflow). These solutions 
do not encounter the magnetosonic curve again. It is also 
possible to construct MHD shock solutions in this scheme. 



crosses the magnetosonic critical curve smoothly and 
therefore we have only one degree of freedom. Numer- 
ical integrations with asymptotic solutions l(4"8"]) and (|4~9"1) 
for the lower signs appear unstable in our experiment. 
This problem remains to be investigated further (Lou & 
Bai 2007 in preparation). 



4.10 MHD Solutions of Inner Thermal 
Fall and Outer Thermal Expansion 

The hydrodynamic thermal fall asymptotic solution at 
small x reported by Lou & Li (2007 in preparation) is in 
a parameter regime quite different from earlier analyses 
(i.e., 7 > 4/3 here). In our polytropic MHD generaliza- 
tion, this asymptotic MHD solution (f3"o| and (|3~7| seems 
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to directly connect with the MHD thermal expansion 
solution ([38 ]) — (|42|) at large x as presented in Figure [H 
While not being a proof, our extensive numerical explo- 
ration seems to indicate that this type of MHD solutions 
does not encounter the magnetosonic critical curve. 

There is one degree of freedom for constructing this 
type of MHD solutions, namely, the choice of inner mass 
point m(0), once parameters 7, n and h have been cho- 
sen. 



5 Discussion 

In this paper, we have shown that more general poly- 
tropic self-similar MHD shock flow solutions can be ob- 
tained under various conditions and they can be applied 
or adapted to various astrophysics systems by examples. 
In particular, we focus on the case of q ^ 0, that is, the 
specific entropy is conserved along streamlines but is not 
necessarily constant in space r and time t. More specifi- 
cally, the specific entropy is related to the enclosed mass 
that varies in both space and time in general. While this 
more general polytropic MHD formalism extends the pa- 
rameter regimes of several known classes of self-similar 
solutions, this same condition of q ^ (or n+j ^ 2) also 
excludes certain solutions, notably the exact global MHD 
solution (Wang & Lou 2007). In astrophysical applica- 
tions, it would be informative and important to know 
these generalizations and constraints for relevant theo- 
retical model development. 

In this paper, we mainly studied quasi-spherical poly- 
tropic MHD gas flows as models for various stellar level 
dynamical processes. However, with various sensible ex- 
tensions and adaptations (such as the two-fluid prob- 
lem with dark matter-baryon matter coupling by gravity 
at the galaxy cluster level; Lou 2005; Lou, Jiang & Jin 
2008; Jiang & Lou 2008 in preparation), similarity so- 
lutions can be applied to different plasma flow systems 
with their own characteristic spatial and temporal scales. 
Another example is the attempt of Suto & Silk (1988) 
to relate self-similar solutions to galaxy formation. In 
principle, if an MHD system experiences a long dynamic 
evolution, it may behave sclf-similarly and thus be de- 
scribed by our model. We emphasize that in Suto & Silk 
(1988), their cases of n = 1 with different values of 7 ^ 1 
do not correspond to specific entropy conservation along 
streamlines but correspond to specific entropy variations 
in time t only. In our cases of n = 1 with different values 
of 7 7^ 1 , the specific entropy is a function of both r and 
t and, in particular, is conserved along streamlines. 

In the preceding section [4j we have systematically 
presented various astrophysical applications and impli- 
cations of our theoretical model development in reference 
to earlier observations and theoretical models in relevant 
astrophysical contexts. These include the following cases. 
(1) The isothermal EWCS (Shu 1977; Shu et al. 1987) 
for the inside-out collapse scenario of star formation can 



now be generalized by our model here in several aspects, 
namely, the inclusion of a completely random magnetic 
field in the collapsing cloud and the more general poly- 
tropic ISM without the restriction of an isothermal gas. 
These two aspects provide the modelling basis for mas- 
sive star formation, radio synchrotron and/or x-ray emis- 
sion diagnostics and radiative synthesis of molecular line 
profiles in star-forming clouds. In particular, a more re- 
alistic non-isothermal temperature profile can be deter- 
mined within a self-consistent polytropic model frame- 
work. This is extremely important for modelling proto- 
stellar cores embedded in molecular clouds. 

(2) The isothermal EECC solutions (Lou & Shen 2004; 
Shen & Lou 2004) for the envelope expansion and core 
collapse scenario of star formation in molecular clouds 
can now be generalized by our model here in several as- 
pects, namely, the inclusion of a completely random mag- 
netic field in the collapsing cloud and the more general 
polytropic ISM without the restriction of an isothermal 
gas. In addition to the implications already noted imme- 
diately above in item (1), we can now take into account 
of inflows (Fatuzzo et al. 2004), outflows of ISM in clouds 
as well as MHD shocks in a more realistic manner. There 
are growing observational evidence for inflows, outflows 
and shocks around star-forming cores in ISM clouds. Our 
MHD shock flow solutions are necessary for certain ob- 
served line profiles (Lou & Gao 2008 in preparation). 

(3) The isothermal magnetized EECC solutions (Yu & 
Lou 2005; Yu et al. 2006) is now replaced by a more 
general polytropic plasma flow which is permeated with 
a completely random magnetic field. By removing the 
isothermal constraint, the radial scaling of mass density 
profiles in a cloud may vary from ~ r _1 to r~ 3 ; ob- 
servationally inferred radial density profiles of molecular 
clouds do fall within this range. For the hot magnetized 
intracluster medium (ICM) in the context of clusters of 
galaxies on Mpc scales (e.g., Lou et al. 2008; Jiang & Lou 
2008 in preparation), we predict the same mass density 
profile range. 

(4) The isothermal shock models (Tsai & Hsu 1996; Shu 
et al. 2002; Shen & Lou 2004) are substantially extended 
by our MHD model here in several aspects, namely, the 
inclusion of a completely random magnetic field, the spe- 
cific entropy conservation along streamlines, and the in- 
clusion of inflows and outflows. Our models are more 
versatile in modelling the so-called "champagne flows" in 
star-forming Hn regions which involve ionization fronts 
and shocks. 

(5) Various hydrodynamic isothermal shocks have been 
extensively explored by Bian & Lou (2005) for possible 
asymptotic flows at large and small radii. For modelling 
"champagne flows" in star-forming Hn regions, Shu et al. 
(2002) considered only central isothermal Larson-Penston 
type solutions. We can now further construct "cham- 
pagne flows" with central free falls, with a completely 
random magnetic field and with MHD shocks. It is also 
possible to construct MHD shock flows associated with 
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different central solutions. In general, rclativistic cosmic- 
ray particles (electrons in particular) can be abundantly 
produced in either MHD wind shocks or accretion shocks. 

(6) The quasi-static model solutions first obtained in a 
conventional polytropic gas (Lou & Wang 2006, 2007; 
Wang & Lou 2007) are now generalized to more general 
polytropic magnetized gas flows in this paper. Based on 
the new quasi-static asymptotic solutions, Lou & Wang 
(2006, 2007) proposed a new class of rebound shocks for 
supernovae. Among other things, we would like to under- 
stand the ejection of stellar materials of the progenitor 
star, the mass of a compact object left behind, as well as 
the origin of intense magnetic fields of compact objects. 
For a rebound shock process in a progenitor star with 
the specific entropy conserved along streamlines, the en- 
tropy distribution in association with the enclosed mass 
becomes an important aspect for various possible solu- 
tions (see also Lou & Cao 2008 for the case of 7 = 4/3). 
This can give rise to a variety of rebound MHD shocks 
in supernova explosions. Of course, we still need to learn 
about the specific energetic corresponding to any chosen 
entropy distribution at a certain stage. 

(7) The asymptotic magneto-accretion solutions with mag- 
netic parameter h > 4 were first studied by Yu & Lou 
(2006) for the case of an isothermal gas and by Wang 
& Lou (2007) for a magnetized conventional polytropic 
gas. This asymptotic MHD solution has now been ex- 
tended to a general polytropic MHD gas flow and may 
be relevant to extremely intense magnetic fields inferred 
for magnetars (e.g., 10 13 — 10 15 G). A strong surface mag- 
netic field of the progenitor star is implied. Magnetic Ap 
stars have surface magnetic fields of several to ten thou- 
sand gauss and give a typical ft- > 4 by estimates. 

While our model results are mainly restricted^ to 
n > 2/3 and 7 < 4/3 for most situations (except that 
the thermal fall and thermal expansion solutions have 
3/2 < 7 < 5/3), the special case of 7 = 4/3 contains a 
variety of substantially new solutions (Lou & Cao 2008) 
that differ from those homologous core collapse solutions 
explored earlier by Goldreich & Weber (1980) in the con- 
text of supernovae. Lou k. Cao (2008) describe this new 
theoretical model development in the absence of a com- 
pletely random magnetic field and obtain various novel 
results as compared to our analyses here. In particular, 
the constant coefficient C in equation (|13p cannot be ab- 
sorbed by a rescaling transformation, and features of hy- 
drodynamic as well as MHD systems will depend upon 
the choice of this constant coefficient C in a nontrivial 
manner. 

Parallel to the analysis of Zhou et al. (1993), in which 
the isothermal model of Shu (1977) was claimed to be 
supported by observational inferences, we are investigat- 
ing radiative diagnostics of the velocity, density and tem- 
perature profiles in collapsing cores of clouds from molec- 
ular spectral line profiles (Lou & Gao 2008 in prepara- 

3 For n > 2/3 and 7 > 4/3, we should proceed with care 
about magnetosonic critical curves with n + 7 > 2. 



tion). Meanwhile, magnetic fields in proto-stellar clouds 
not only affect their dynamic evolution, but also pro- 
vide diagnostic signals, such as synchrotron radiation 
and shock acceleration of relativistic electrons etc. By 
including radiative transfer processes, our model is ca- 
pable of predicting such diagnostic signals. 

We also note earlier theoretical analyses on the insta- 
bility of known self-similar solutions (e.g., Ori & Piran 
1988; Hanawa & Nakayama 1997; Hanawa & Matsumoto 
1999, 2000; Semelin, Sanchez & de Vega 2001; Lou & Bai 
2008 in preparation). Instability analysis helps to test 
whether certain self-similar flows will actually occur in 
nature; or if they do occur, whether they will last long 
enough. Our model provides a variety of new polytropic 
MHD self-similar solutions with or without shocks, all 
ready to be examined for such an instability analysis. At 
this stage, these stability questions remain completely 
open. 

Finally, numerical simulations are powerful and nec- 
essary in determining how likely our various self-similar 
MHD solutions will eventually emerge in a sensible way 
for plausible initial and boundary conditions based upon 
astrophysical input. Moreover, nonlinear instability anal- 
ysis also needs to employ numerical simulations. We hope 
that the results presented here will trigger extensive nu- 
merical explorations for transient behaviours leading to 
these self-similar MHD solutions with specific entropy 
conserved along streamlines. 
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A Comments on the Positive Definiteness of 

nx — v 

For a self-similar gas flow where the self-gravity may or may 
not play a significant role and for the self-similar transforma- 
tion in the form of 



.1/2 



r — k t n x , u = k t 



l/2,n-l 



p - a/{AiiDt s 



(52) 



where D and k are two constant positive parameters and n 
and s are two exponents with 3n — s > 0, we then derive from 
mass conservation equation ([lj 



M(r 2 ) - M(ri) = 7- rR ax ( nx ~ v ) 

(6n — s)JJ 



(53) 



where Xi = r i /(k 1 ' 2 t n ) with subs crip t i — 1, 2 at a given 
time t. The condition for equation (|53p to be valid is that the 
gas flow within the radial range between ri and r\ evolves in 
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Fig. 6 Semi-complet e M HD solutions with inner thermal fall 
asymptotic solutions (|36|l and (|37p at small x and outer MHD 
thermal expansion asymptotic solutions at large x (see sec- 
tion 3.1.5) for 7 = 1.65, n = 0.8, q = 2.25 and h = 0.5. The 
'LL1' and 'LL2' MHD solutions have parameter m(0) — 0.4 
and m(0) = 0.2, respectively, and are constructed by inte- 
grating from small x outwards with above parameters. The 
outer portions of these two MHD solutions converge to the 
same asymptotic solution at large x to the first order, with 
the proportional constant c of v(x) ~ ex being 0.5063 for the 
thermal expansion solution and the corresponding constant 
coefficient E for this solution is 2.277 ', and the value of P as 
the power law index for a(x) versus x is —1.638. By numerical 
experiment, these two different solutions at small x gradually 
merge to the same asymptotic solution at large x and do not 
encounter the magnetosonic curve. 



In contrast to the above case, we now consider the sit- 
uation where the self-gravity is indeed important. If we are 
to establish a self-similar evolution, the total enclosed mass 
M(r, t) would enter the radial momentum equation and thus 
be a self-similar variable under consideration. From a dimen- 
sional analysis, if we take only k and D as dimensional pa- 
rameters, the sensible transformation for total enclosed mass 
M(r, t) in the power law form should be defined as 



M 



u3/2^3n-s 

'/ 

(3n - s)D 



i(x) , 



(55) 



where the reduce d en closed m(x) depends on x only. Com- 
paring equations (|54|l and (155[) . we find (reducing the overall 
scaling factor) 

m(x) — m(xo) — ax (ni- v)\% , (56) 

or we have 

m(x) = ax (nx — v) + constant. (57) 

However, since equation ([T]) only takes care of regions where 
r > while the mass at the origin r = is not included, 
this constant is not arbitrary. A difference in this constant 
means a difference in the total enclosed mass that changes 
with time in a scaling of t in ~ s . The equation that includes 
the mass at r — is equation ((2)1 , and this fixes the constant 
in equation (|57|l to be zero. Thus, in order to introduce the 
total enclosed mass into the radial momentum equation as a 
self-similar variable, nx — v > is required for 3n — s > 0. 
Otherwise for 3n — s < 0, we should require nx — v < 
accordingly. 



B Determination of the Magnetosonic Critical 
Curve 

The magnetosonic critical curve of equation (|15p is defined as 
the curve along which both the numerators and denominators 
{ A(x, a, v), X(x, a, v) and V(x, a, v) } of equation (|15[) van- 
ish. In fact, only two of these three vanishing conditions are 
independent and we shall simply use the intersection of the 
two surfaces A(x,a,v) = and X(x,a,v) — to determine 
the magnetosonic critical curve. 

In order to obtain the magnetosonic critical curve nu- 
merically, we denote 9 = nx — v and 9 > is required for 
a positive total enclosed mass (see Appendix [A}. We readily 
obtain from A(x,a,v) = X(x,a,v) — the following two 
equations, namely 



a self-similar manner. Equivalently, we can express the total 
enclosed mass M(r) as 



, rfl . I 2 - n a 
a= ^+( 3^2 + 2 



(n — l)(nx + i 



■20'' 



h.3/2_^3n — 8 

M{r)= (3n-s)D ax2{nX - v) 



r/ik 1 / 2 ^) 



7- /(fcl/2t") 



+ M(ro) , (54) 



as long as the gas flow evolves in a self-similar manner within 
the radial range between r and ro. 

In general, nx — v > is not necessarily required be- 
cause M(ro) may not behave self-similar ly (in tim e or in 
space); thus the only requirement for equation (|54[) is that 
ax 2 (nx~ v) is a monotonically increasing function of x, which 
is guaranteed by the mass conservation equation, provided 
that a > and x > 0. This is why in self-similar flow sys- 
tems where the self-gravity is not important (e.g., Chevalier 
1982), solutions with nx — v < are allowed (e.g., for the 
reverse shock portion of the solutions in Chevalier 1982). 



and 



qhx 



2 - U < q ^(2hx9- 



3n - 2 + 2 



3n-2 



(58) 



2 _ ,, + ^l q ) a i-"+3" 8 /V 9 8 + hax 2 = 9 2 . (59) 



A direct substitution of equation (158[) into equation (159[) to re- 
place a yields a nonlinear equation for x and 9 which leads to 
a curve 8 = 6(x) for the magnetosonic critic al cu rve. Using 
this curve, the definition of 9 and equation (|58p . the mag- 
netosonic critical curve is then determined in the form of 
(x,a,v). In our analysis, the relation 9 = 9{x) is determined 
numerically for a set of three specified parameters {n, q, h}; 
the definition of q contains the polytropic index 7 parameter. 
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C Eigensolutions crossing the Magnetosonic 
Critical Curve 

Using the l'Hospital rule and with a superscript prime '/' 
indicating a differentiation with respect to x, we obtain the 
following differential relation along the magnetosonic critical 
curve for v' and a' from equation X(x,a,v)v' = V{x, a, v), 
namely 

[fi(x,a,v)a' + f 2 (x,a,v)v' + f 3 (x,a,v)]v' 



+ 



. . 2(x-v)a t , . 

xinx — v) 

f 4 (x,a,v) - f 5 (x,a,v) 



(nx — v) 
2(x — v)a 



+ 



x(nx — v) 



fi(x,a,v) - f 6 (x,a,v) =0 



(63) 



The two real roots v' of quadratic equation (|63p on the mag- 
netosonic critical curve represents two eigensolutions which 
can go across the magnetosonic critical curve smoothly. 



= f4,(x, a, v)a + fs(x,a,v)v' + f 6 (x,a,v) 



(60) 



where functional coefficients fi(x, a, v) with i — 1, • • • , 6 are 
defined by 



3nq\ 



f-L(p,a,v) = hx 2 + ( 2-nH — q\\l-n+ y ] 



xa- n+3nq/2 x 2q 8 q , 
f 2 (x, a,v) = 2e-q(2-n+ ^-^q 



xa l-n + 3n q /2 x 2 qeq -l ^ 

fs{x,a,v) = 2hax - 2n6 + 2q( 2 -n+ -^— — q 



., l-n+3nq/2 2q-l nq . In, 3n — 2 

xa H! x H 6 H + nql 2 — n -\ q 



xa 1 - n+3nq/2 x 2q e q - 1 



fi(x, a, v) = 2(1 — n)hx 



(3n - 2) 

+ [2(1 - n)(2 - n) - n(3n - 2)q](l - n + 3nq/2) 
xa~ n+3nq/2 x 2q $ q + [2(2 - n) + (3n - 2)q] 

x(l - n + 3nq/2)a- n+3nq/2 x 2q - 1 e q+1 , 

2a8 



f 5 (x,a,v) = (n - l)(nx - 28) + 



(3n - 2) 



-q[2(l-n)(2-n)-n(3n-2)g]. 
-(g+l)[2(2-n)+(3n-2) g ]a 



i-^nVinq/2 2q nq-l 



1— n43n<j/2 2q-l nq 



fe(x, a, v) = 4(1 — n)hax — 2n (n — l)x + n(n — l)v 
+2n(n - 1)6» - 2na6/(3n - 2) 
+2q[2{l-n)(2-n)-n(3n-2)q]a 1 ^ +3nqft x 2q - 1 6 q 
+nq[2{l-n)(2-n)-n(3n-2)q]a 1 - ni - 3nq/2 x 2q 6 q - 1 
+ (2q-l)[2{2-n) + {3n-2)q]a 1 - n+3 " q/2 x 2q - 2 6 q + 1 



D Solution to MHD Shock Conditions 

In reference to MHD shock jump equations (|22|1 — (1241) . we 
may denotsQ A = n — Vi/xt and rearrange these MHD shock 
equations into the following equivalent form 



aiAi = 02/2 



2-n+3nq/2 3q-2 r q 



•l 



x?-*I? + <*i A + 



2 — rc.+3nq/2 3<j — 2 r-iq 



ha 1 



x a 2 q -r q + aa r; + - 

27 ^ n+3ng/ 2 x 3 q -2 r , + r 2 + ^^ 

(7- 1) 

_ _27_ l-n + 3nq/2 3q-2 r q , p 2 



(7-1) 



x% q ~'r q + A + 2ha2 



(64) 



Using the first mass conservation equation in (|64f) , the down- 
stream reduced density 01,2 can be immediately replaced to 
obtain 

2-n+3nq/2 3q-2 r q . n 2 hai 



'a^ A + «i A + 



2 — n+3n<j/2 p2 — n + 3n<j/2 



A 



Xo 9 +aiAA + — 

2 — n+(3n — 2)q/2 * - > 



ft a? A 2 



A 

27 a i_„ + 3„,/2 a .3,-2 r , + r 2 + 2hai 

(7- 1) 

,-. 1 — n+3ng/2 t-,1 — n+3nq/2 

_ ^7 Q l A 3q-2 

fry _ 1) pl-n+(3n-2)«/2 2 



2 A 2 ' 



+A + 2ft 



aiA 
A 



(65) 



Eliminating x 2 q ~ term in both equations of (|65[) . we obtain 
a cubic equation for A of the downstream side (subscript 2) 
in terms of upstream variables (subscript 1) 



+n(g+l)[2(2-n) + (3n-2)g]a 1 ^ l+3n ' l/2 a; 2 ^ 1 6» 9 



(61) 



[see equations (|15p and (|16p ]. Together with the differential 
relation 



, av' — 2(x — v)a/x 

a = — ■ — 

(nx — v) 



(7 + l)i-i3 ( l-n+3nq/2 3<j-2 n q-l . n . ft«l \ , ,2 
-1 2 — I CKi Si i, +^H J 2 

27 V 2A ' " 



l_„+3n 9 /2 3 9 -2 r q , 7~ 1 r 2 , 7_~ 1 



(62) 



•a.-"-^ + -!—^n + -! ^ftai T 2 

2 7 7 



we obtain an algebraic quadratic equation for v 1 = dv/dx at 
the magnetosonic critical curve, namely 



27 



(66) 



(nx — v) 



fi(x,a,v) + f 2 (x,a,v) 



(v'f 



4 Please note that in defining similar new variables above 
equation (46) of Lou & Wang (2006), there was a typo in this 



regard. 
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Finally, b ecau se A = A is a trivial and unphysical solution of 
equation (|66l) for an MHD shock, we then obtain a quadratic 
equation of A for possible physical MHD shock solutions 



(7+l) n2 

o •* 2 

2 7 



1— n+3nq/2 3q-2 pq-1 



rr 1 + 



7' 



27 



-A 



2A 



7) 



27 



/iQl = , 



A 



(67) 



and the positive A root of this quadratic equation represents 
a proper MHD shock solution. 

For A > and 1 < 7 < 2, quadratic equation (I67|l al- 
ways has two real roots for A with the plus-sign root being 
positive and the min us-s ign root being negative. The first 
relation of equation f|64f) gives the corresponding 02- One 
can usually determine a sensible X2 and thus V2 and A val- 
ues. Occasionally, we encounter the unrealizable situation of 
a complex 22 caused by the presence of a random magnetic 
field (i.e., /i / 0). For h — 0, the situation of a complex X2 
would not arise. This suggests that self-similar MHD shocks 
cannot form within a certain radius and time for a given set 
of parameters. 
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